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Abstract: In order to prepare the ground for evaluating classes of three-loop sum-integrals 
that are presently needed for thermodynamic observables, we take a fresh and systematic look 
on the few known cases, and review their evaluation in a unified way using coherent notation. 
We do this for three important cases of massless bosonic three-loop vacuum sum-integrals 
that have been frequently used in the literature, and aim for a streamlined exposition as 
compared to the original evaluations. In passing, we speculate on options for generalization 
of the computational techniques that have been employed. 
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Introduction 





Our knowledge about sum-integrals, needed for evaluating phenomenologically relevant equi- 
librium observables in thermal field theories (some examples being [1-3]), is by far not as 
developed as the knowledge about continuum integrals [4] needed for standard high-energy 
(but zero-temperature) phenomenology, as e.g. reviewed recently in [5]. One of the reasons 
seems to be the impenetrable structure of the multiple infinite sums that are involved. 

Even considering the simplest class of sum- integrals, dimensionally regularized massless 
bosonic vacuum sum-integrals (which constitute zero-scale problems since their only dimen- 
sionful scale - the temperature T - scales out trivially) that we will call hot tadpoles in the 
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following, only a very limited number of cases are known in practice (for a review, see [6]), let 
alone the number of technical tools that have been developed for handling such sum- integrals. 

All such presently known hot tadpoles contain one-loop two-point sub- integrals n(P), 
whose structure is heavily exploited in the process of evaluation. While the 1-loop tadpoles 
can be computed exactly, i.e. as functions of the dimension d, all 2-loop sum-integrals can 
be factorized into 1-loop ones e.g. by systematic use of integration-by-parts (IBP) methods 
[5, 7], the first non-trivial cases occur at the 3-loop level. 

It turns out that the few hot 3-loop tadpoles that have been evaluated in the literature 
have been dealt with on a case-by-case basis by hand, often in some painstaking process, 
involving inspired tensor transformations, elaborate (UV- and IR-) subtractions, skillful inte- 
gration tricks, mixed momentum and coordinate space techniques, numerical integration and 
the such [8-10]. 

At present, however, there are not only pressing open three-loop questions [11], but even 
interesting open problems at four-loop order [12] that involve a (large) number of yet unknown 
sum- integrals. While exactly one type of 4-loop tadpole has been evaluated so far [9], it is 
clear that a more systematic treatment is urgently needed. The purpose of the present note is 
to re-analyze the known non-trivial 3-loop cases and to streamline their derivations in terms 
of a unified notation [13], in order to prepare the ground for tackling further 3-loop and 4-loop 
tadpoles in an efficient way. 

The relevance of hot tadpoles can be appreciated from the following: In modern treat- 
ments of equilibrium thermodynamics, despite the known problem of infrared (IR) divergences 
in massless thermal gauge theories [14], an effective field theory (EFT) approach [15] allows 
for clean separation of IR and ultraviolet scales; the sum-integrals treated here constitute, 
in the jargon of those EFT's, the hard contributions to the corresponding observables, and 
hence need no IR regulator (other than being dimensionally regularized) [16]. 

The three key examples of known non-trivial massless 3-loop vacuum sum-integrals that 
exhibit a range of useful techniques are of 'spectacles-', 'basketball-' and 'tensor-spectacles- 
type', and can be represented as special instances of a general class as defined by (cf. 
Eqs. (A.23),(A.30) of Ref. [2]), 

Our Euclidean notation is such that we use bosonic four-momenta P with = pg -|- = 
(27rnpT)^ -|- p^, and where the sum-integral symbol is a shorthand for 

with d = 3 — 2e and the sum taken over all integers Up G Z. In the notation of Eq. (1.1), A1o,o 
and A^i,o stre the basketball- and spectacles-type 3-loop tadpoles, respectively, while we refer 
to A^2,-2 as the tensor-spectacles-case, owing to its numerator structure. 
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In the remainder of the paper, we wUl discuss the cases ^Alfi, A4o,o and A^2,-2 of Eq. (1.1) 
in turn, in Sections 2-4. Sec. 5 contains conclusions, while a few technical details are relegated 
to the appendices. 

Before turning to the specific cases, let us note that (by splitting AQR = {Q + Q){R + R) 
and exploiting the denominator's invariance by shifting Q —Q — P or — )• —R — P in the 
second instance only), the sub-class 



is seen to involve scalar 1-loop sub- integrals 11 only . It can hence be treated with the scalar 
methods employed for ^Alfi and Alo,o see Sec. 2 and Sec. 3, respectively. As an example, 
from Eq. (1.3) we immediately get Aii^-i = — 7Wo,o/2, where we have used that the 2-loop 
sunset sum-integral 



vanishes identically in dimensional regularization as can be shown via integration-by-parts 
(IBP) techniques [5, 6]. 

2 The 3-loop spectacles 

In the notation of Eq. (1.1), the 3-loop spectacles-type sum-integral A4i.o is defined in terms 
of the 1-loop 2-point function 11 as 



We shall translate the original computation of A^i,o from Ref. [10] (relying on the methods 
pioneered by [8]) to our notation and systematics, using less than their seven pages, in a 
transparent way, and in a notation that is generalizable to other cases. 

2.1 Decomposition of A4ifi 

The spectacles can be identically re-written as 




(1.3) 



(1.4) 




(1.5) 




(2.1) 



+ 




To unclutter the text, we have coUected various definitions in the Appendix (cf. Eq. (A. 2) for Ji). 
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where Sp^ picks out the Matsubara zero-mode, the primed sum excludes the zero-mode and 
we have suppressed the argument {P) of ah functions in curly brackets. This re-written form 
becomes useful if TIa^b,c have the structure ^ifi{T,e) / and Hd = f{T,e) does 

not depend on the momentum P. For then, the first term of Eq. (2.2) is proportional to the 
2-loop sunset sum-integral 5, which vanishes as already discussed above; all terms that do not 
involve 11 are trivial 1-loop tadpoles /, which are known analytically in d = 3 — 2e dimensions 
(see Eq. (A. 2)); the zero-mode term involving [Ilyi — Hd] n is known analytically in terms of 
the 2-loop function A introduced in Eq. (21) of [13] (see Eq. (A. 3)); and in the three remaining 
terms that involve 11, the three subtraction terms ^a.b.c can be independently chosen such 
as to facilitate their evaluation. 

The strategy of Ref. [10] amounts to choosing^ (for a motivation of this choice, see below) 

/3 rG(l,l,rf) p /3 /3 

where j3 = G(l, l^d + 1) with G given in Eq. (A.l) and a is a constant to be fixed later. 
Eq. (2.2) then reduces to 

TWio = + ^ + 2/3^(1+^,1,1) + 2rG(l,l,(i)A(3/2 + e,l,l) -0 + 

+ B + 2C- P^h+2e + 2p^I^ + APhI-^ , (2.4) 

where we have introduced the shorthand f{x) = f{x + e) — f{x)/{a T'^Y for convenience. The 
first zero in Eq. (2.4) is the 2-loop sunset discussed above and the second zero comes from 
the fact that 2^ Sp^ \^a\^ /P'^ is scale- free and hence vanishes in dimensional regularization. 
We have introduced the notation B, C for the three non-trivial 3-loop sum-integrals^ 

A^^J^[U-Ua]\ B^^^^[U-Ub]\ C^^^^[Ub-Ud][U-Uc] (2.5) 

that involve IT, and whose evaluation we shall discuss in the next subsection. In fact, the 
subtraction terms TiA,B,c,D defined above were chosen such that A,B,C are finite in d = 3, 
and can be evaluated numerically after simplification by e.g. the spatial Fourier transform 
method of [8]. More concretely. Ha subtracts the leading UV- and IR-divergences in A; Hb 
subtracts the leading UV-divergence in ;B; n^i was chosen such that [Hs — n/5] is finite as 
d — > 3; and lie subtracts the leading and sub-leading UV-divergences in C. 

2.2 Evaluation of A, B, C 

Let us now bring the 3-loop sum-integrals A, B, C into a form suitable for numerical evaluation. 
The (inverse) 3d spatial Fourier transforms that will be used below read [9] 

7^} = (If Z^-'" I"'" p-^-' - 2^'^} ^ 

^We take a slightly different Tin here, whose leading term at e — >■ equals the choice l/(167r^e) of [10]. 
^In the notation of Ref. [10], A = J", B = If and C|^=i6^2/e7E = ^3 l^=4^T/e(i+^E/2) ■ 
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For A, re-writing 5^, [H - Ua] = Sp, [U-Ub]-1^x + ^ie) (where i = 

G(l,l,3), while the extra integral is unity and introduced here for notational simplicity); 
using the 3d spatial Fourier transform of [11 — H^] (at pQ = 0); integrating over angles via 
i J\ d-ue^P™ = ^i^^; and letting |r| = x/{2ttT), \r'\ = y/{27rT), |p| = 2ttTz: 



(47r)4 
°° dx 



A + 0{e) , 



A 



A = I — ( coth(x) - - - 1 
— coth(x) 1 

X 



(2.7) 



X 

~ dx 
X 
~ dx 

X 

~ dx 





— [ coth(y) 
y 

— ( coth(y) 
y 



coth(x) 
coth(x) 



1 

X 

1 

X 



dy 
y 

sinh(x 



coth(y) 



y / TT 
1 

- - 1 

y 

1 







°° dz 

sin(zx) smizy) 



\x + y\ - \x-y\) 



1 (2y) 



In 



In 



sinh(x) 



9.5763057898... . 



For jB, using the Fourier transform of [H — 11^]; recognizing /pe^P"" 



.-IpoI'" 

^4^; sum- 



ming over pq via geometric series; scaling |r| = x/{2ttT), \r'\ = y/(27rT); and integrating over 
angles: 



B 



dx 







1 



dy 



— ( coth(x) ) / — ( coth(y) 



1 



In 



sinh(x -|- y) 
sinh(x) 



(2.9) 

(2.10) 
(2.11) 

For C, using the Fourier transform of [11 — He]; expanding [He — H^)] = In -|- 
0{e); integrating over angles; letting |p| = \pQ\y, |r| = x/{2irT); and using the exponential- 
integral Ei(x) = - ^ f f e'^' /o" 1^ TO = e''^lEi(-2|z|) + 7e + In ^: 



0.058739245719... . 



1 



3 (47r) 



C + 0(e), 



c.-of^(».M.)-i-|)| 



Ei(-2nx) + e 



-2nx 



In 



2x ae^E 
n IGvr^ 



0.003496... . 



(2.12) 

(2.13) 
(2.14) 



where the numerical value is given for a = IGir'^/e'^^ and corresponds to^ Eq. (D.27) of [10]. 
For an discussion of the numerical evaluation, we refer to App. B. 



^Note that for this choice of a, we avoid the computation of ^ in Eqs. (D. 20-25) in [10]. 
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2.3 Result 

Expanding Eq. (2.4) around d = 3 — 2e (for a = IGvr^/e'^^), we finally obtain 

Mi,o = ^ ^ ^^-^ [l + vie + V2e^ + ^(e^)] , (2.15) 

4 C'f-1) 

^;i = 3+47E + 2^^^, (2.16) 



1 

V2 = — 

3 



46 - 167E^ + + 241n2(27r) - IO471 - 87E - 247e ln(27r)+ 



c'(-i) , o.c'(-i) . oC"(-i) 



+ 167ETf-rT + 24^rf-f + 2 



38.5309... , (2.17) 



C(-i) C(-i) C(-i) 

which coincides with Eq. (D.51) of [10]. The numerical value in Eq. (2.17) is —4 (yl + S — | C) . 
3 The 3-loop basketball 

In the notation of Eq. (1.1), the basketball-type sum-integral A^o,o is defined in terms of the 
1-loop 2-point function 11 as 



^A^o,o^^[n(P)]2 . (3.1) 



Historically, the evaluation of A^o,o was performed by in Ref. [8], where many of the techniques 
that were later generalized to other cases of sum-integrals, such as basketball-type tadpoles 
with different powers on the propagators and/or factors in the numerator [11, 13], were 
introduced. Here, we translate this pioneering computation of A^o,o to our notation and 
systematics. 

3.1 Decomposition of 7Wo,o 

The basketball can be identically re-written as 

-^0,0 = 1^ {2n,, n} + ^ 5p„ {[n-nB]2 + 2 [Hb-Hb] h - [Hb]"} + ^'{2 [Hc-Hb] n} + 

+ ^'{[n - He]' + 2 [Hs - Hd] [H - He] - [Hcf + 2 [Hs - H^] H^} , (3.2) 

where, in full analogy to Eq. (2.2), picks out the Matsubara zero-mode, the primed sum 
excludes the zero-mode and we have suppressed the argument (P) of all functions in curly 
brackets. 

Let us slightly refine the strategy of Ref. [8] by choosing 

nB=(^. nc = HB+-^, i[o = j-^^^ nB = nc+8|p2|2 — jpi — ' (3.3) 
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where /? = G(l, 1, d+ 1) as above, a is a constant to be fixed later and J„ = ^4^3"/2r(^^y^"^ 
as in Eq. (B5) of [8]. With this choice, the first term of Eq. (3.2) can be shifted to the square 
of a trivial 1-loop tadpole /; for the last term in the first line of Eq. (3.2), note that after 
Y^' — > X] — X] 5po f^^l snm is proportional to the 2-loop sunset sum-integral S and hence 
vanishes; the rest as well as the other zero modes involving [n^ — II/)] 11 and [lis — He] H 
are known analytically in terms of the 2-loop function A of Eq. (A. 3); all terms that do not 
involve 11 are trivial 1-loop tadpoles /; and in the three remaining terms that involve 11, the 
subtraction terms IIb,c,d,e have been chosen such as to subtract UV divergences in order to 
render the sum- integrals finite. 
Eq. (3.2) then reduces to 

Mo,o = 2/3/i/i/(ar2)^ + P + 2/3^(0, 1,1) -0 + 0-4/1^(1,1,1)+ (3.4) 
+ f + 27- - - 4/3/1 h+e -4/1/1/2 + 2/32 /, + 4/3/1 /i + [/2-(a!+l)/|] , 

where we have again used the shorthand f{x) = f{x + e) — f {x) / {aT'^Y . The first zero in 
Eq. (3.4) comes from the fact that ^ Sp^ Pb]^ is scale-free and hence vanishes in dimensional 
regularization and the second zero is the 2-loop sunset discussed above. We have introduced 
the notation T),£,T for the three non-trivial 3-loop sum-integrals^ 

v = ^^5p,[u- Ub]' , ^ ^ ^'[n - ^c? , J- = ^'[n^ - Ud] [U - Ue] (3.5) 

that involve 11, and whose evaluation we shall discuss in the next subsection. 
3.2 Evaluation of V, £, T 

The 3-loop sum-integrals P, T are finite as d 3, such that the 3d spatial Fourier transform 
method of [8] proves fruitful. The (inverse) transforms needed below are [8, 9] 

{n, Hb, He - Hb, - He, - nz)} = 

^ ^d'%^pre-M.L,th(r-) + M, M + i, i±M:Uo(,), (3.6) 



(47r)2 y r2 I 27rr' ItxT f' 3' 45' 

where f = lirTr, and for the last term we have expanded [IIb — n/)] = In + 0{e) 
and used the inverse Fourier transform of the logarithm as derived e.g. in Eqs. (D.11),(D.12) 
of [8] 6. 

For D, using the 3d spatial Fourier transform of [11 — IIb] (at po = 0); integrating over 
angles; and letting |r| = x/{2ttT): 



(4^)2 



T^ = 7JZy2D + 0{e), (3.7) 



'^In the notation of Ref. [8], O = Eq. (2.34); £ ~ Eq. (2.31, 32); T ^ la as treated in Eqs. (D9-D14) of [8]. 
®Note that it is defined up to a Delta function 5{v) which however vanishes in T below. 
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D 



( cotn(xj 







x^ 



2C(3) 



(3.8) 



where the analytic value was obtained via the recursion of App. C. 

For £, using the Fourier transform of [IT — He]; integrating over angles; summing over po 
via geometric series and re-writing 2/(e^^ — 1) = coth(x) — 1; and scaling |r| = x/(27rT): 



£ 

E 



(47r)2 

1 dx 

Jo X 

1 



E + 0{e), 



2 L ^2 \ coth(x) 



(coth(3;) 



1) 



18 



67 , , , 36C(3) C'(-3) ^C'(-l) 



(3.9) 
(3.10) 
(3.11) 



where the analytic value was again obtained via the recursion of App. C. 

For J-", using the Fourier transforms of [11 — II^;] and [11^ — Hd]', integrating over angles; 
letting pq = 27rTn, |r| = x/{2ttT)] and summing over n via geometric series and re-writing 
2/(e2^ - 1) = coth(x) - 1: 



F - 



F + 0{e) 



^ I dx 
F = — 



1 X x^\ 



X 



1 

180 



23 90C(3) a-3) O-l) 

--7- + 37E + — + 2 - 5 

6 TT^ C(-3) C(-l) 



(3.12) 
(3.13) 
(3.14) 



where the analytic value was obtained by first introducing the regulator x^ as in App. C; re- 
writing 2 coth(x) coth(x) = (Sa; coth^ (x) ; integrating by parts all terms involving dx while 
dropping all boundary terms, which vanish due to the regulator; and using the recursion of 
App. C, letting 5 — >■ in the end. 



3.3 Result 

Expanding Eq. (3.4) around d = 3 — 2e {a does not contribute yet), we finally obtain 
^4 (47reTET2)"^' 



Mo,o 
hi 



(47r)2 24 e 
37 327E 



[l + hie + Oie^ 



9 15 

_ 91 C^(-3) C'i-1) 

15 C(-3) C(-l) 

which coincides with Eq. (2.36) of [8]. 



+ 81n(27r) - + — ^^^^ + 24(D + E + 2F) 



7r2 15 C(-3) 



(3.15) 



(3.16) 



- 8 - 



4 The 3-loop tensor spectacles 



A first non-trivial representative of the class Eq. (1.1) involving numerator structure is A^2,-2- 
In this case, we do not have an easy way of dealing with the scalar products in the numerator 
as was still the case for Mi-i, see Eq. (1.3). Here, we wish to first relate the computation of 
A^2,-2 to an auxiliary one [2], and then evaluate the latter, for historical reasons denoted Xgqed 
[8], using our notation and systematics. Let us note that ^^2,-2 is quite another category 
compared the previous two, as it needs tensor methods. Let us re-write'' 

4^2,-2 = 4(5 - d)/2/i/i + Mo,o + 3^ [^^..{P)]^ (4.1) 



p 



with U,^{P) = 2hg,^-f^ QHQTW ' 

where I[^u{P) was chosen transverse, P^n^jy(P) = 0, consequences of which will be exploited 
next. 

4.1 Relating M2-2 to Zsqcd 

Transversality P^Il^i^{P) = constrains the structure of the symmetric tensor 11^,^ to 

n^,(P) = A^,Ua{P) + B^,Ub{P) , (4.3) 

where A^u = Au^, = are projectors A A = A, BE = B, AB = (with traces 
tr^ = d — 1, tr B = 1) which are orthogonal to the external momentum PA = = PB. 
Concretely, trading the 4-vector U = (1,0) for the linear combination V = P^U — {PU)P 
that satisfies PV = 0, we have A = g — -pr — (for which in fact also VA = = UA) and 
B = X^. The scalar coefficients fl^x^B of Eq. (4.3) can hence be obtained via projections^ 



V 

3212 



^ ^ ^ ^ (- - ^-"-0 ^ ^^^^^^^ ■ 



tr^n(P) 


1 


tr^^ 


d-1 


trBU{P) 




tvBB 


y2 



^"^^^ = ^^trm ' = f^Mf^-n^-(^) = ^ noo(P) . (4.5) 

For the specific case at hand, we can read off H^^ and fioo from Eq. (4.2) to obtain 

n,, = p2n(p) + 2/i(d-i), n^P) = ^(2/1 - q?^qI%^ ■ (4-6) 

Returning to Eq. (4.1), A^o.o is the 3-loop basketball sum- integral, while the last term is 
related to 2sqed of Appendix H in [8] by an IR subtraction in the po = mode defined by 

t^^^^A^J^l^B^J^^ (4.7) 



^Note that in our conventions d = 3 — 2e and hence the trace of the metric tensor is g^^^ = rf + 1. 
^Clearly, Hg will be the "hard" case whenever it occurs, involving 1/p^ and go etc. 
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when choosing^ IV^ 



9 rrlR — 2(d 2)li+4l 2_ -Qffi, _ 2j_^ _ momentum-independent: 



-^sqcd 



n„„ + 2 fn 



IR 



2 MM 



Hb ~I~ Ogcale— free 



Xsqcd + 2nif 1, 1; 0) + 8 (njf - n^B^) ^(2, l, l; 2) + o . 



(4i 



Here, the 1st hne is a trivial re- writing; the 2nd Hne follows via Eqs. (4. 3), (4. 7), using the prop- 
erties of the projectors, plugging in Eq. (4.4), noting that n|^g| are momentum-independent 
and dropping scale-free integrals that vanish in dimensional regularization; and in the last 
line we have used Eq. (4.6), again dropped scale-free integrals and expressed the remaining 
sum- integrals in terms of the 2-loop tadpole A from Eq. (A. 4). 

4.2 Evaluation of Xsqcd 

Owing to Eqs. (4.1) and (4.8), instead of 7^2,-2 the authors of Ref. [8] choose to compute 



-^sqcd — 



5 n^^l 



uv 



5 n^^l ^ + flu^ r2n 



} 



(4.9) 
(4.10) 



d-l 



d-1 



+ 



n 



uv 



+ jp2 ( ^n^M - ^'nuv _ 26,, [id - m'^ + 



IRl 



dUB -2Ua 



+ 



rUV 



2U 



p2 



n 



uv 



+ 0, (4.11) 



where in the 2nd line a UV subtraction nJ/J = {A^.^ + B^„) ^n^v 

was introduced^*^; for the 
3rd line we have used projector properties as well as Eq. (4.4), and defined 

Ha = ^ (n,, - p^nuv _ 5,, [id - + ) , 



p2 

1 
p2 



(4.12) 
(4.13) 



and for the 4th line we have assumed = Ci/[P^]"' and dropped scale- free integrals. 

The essence of the computation of Ref. [8] is now the treatment of the terms involving ILb 
in Eq. (4.11) (note the similarity to Eq. (H.14) of [8]). Choosing different UV subtractions 
for the zero- and non-zero modes via 

= Ub + il - Sp,)d2Ii/P^ these terms are finite 
and can be treated in d = 3 by spatial Fourier transform methods. Ref. [8] states the simple 
result 

1 /■ - r - - 1 /■ i2 



d-2 



P 



dUB -2nA 



Ha 



+ Oie), 



(4.14) 



^This particular choice in fact reflects IIi^ gj — II{a,b}{0, p — >■ 0). 

^°Note that its tensor structure is {A^i, + B^v) ~ (ff/ji/ — P^Pu/P^), as expected at zero temperature. 



- 10 - 



which follows from a rather lengthy calculation, involves an "amazing cancellation" and is 
explained in our Appendix D. As a result, the computation of Xsq^j (up to the constant term) 
is reduced to elements that already appear in the basketball case A^o,o- 
In detail, with Eq. (4.6) and the choices of fl^^ and given above, 



UA = U-UB-{l-5p,)2h/P^ , 

nB = ^ (iib - p^ns/d - 2/1 + 6p,4il) . 



Plugging Eq. (4.14) into Eq. (4.11) then results in 



^qed 



Ha 



\-,T li- 



ny I 2n^^ 



p2 



n 



uv 



0{e) 



(4.15) 
(4.16) 



(4.17) 



= ^ {[n]' + ^([4/i/p2+n5_2n] + {l-6pJd[2h/P^]')] + Oie) . (4.18) 
Recognizing in Eq. (4.18) the term quadratic in 11 as Mq^, re-writing the linear term as 

^ {nns} = [n - n^,] [n^ - n^] + n^, [n^ - n^] }+ 



(4.19) 



where the first term on the right-hand side (rhs) is J- (of the A4o,o calculation of Sec. 3, cf. 
Eqs. (3. 5), (3. 14)) and the others are elementary (as are the remaining terms in Eq. (4.18)), 
2^sqed evaluates to (using again /(x) = f (x + e) — f (x) / {aT'^Y as well as the functions collected 
in Eq. (3.3) and in App. A) 



d-1 



Isdcd = Mo,o + —- /2e + 4/3/1 h+e - 2\ T + (3' I, + 2 ^ h + 



8/3 Ji 



+ /3 ^(0, 1, 1; 0) + /3 Ii h/{a T^)^} +Adhhh+ Oscale-free j + 0{e 

The expansion around (i = 3 — 2e coincides with Eq. (H.30) of [8]: 
(47reTET2)"^' 



[/2-(d+l)/|] + 

(4.20) 



^sqcd 



108 



23 

2^ + 



+ 0{e) 



(47r)2 
4.3 Result 

Putting together Eqs. (4.1), (4.8), (4.20) and expanding around d = 3 — 2e, we finally obtain 



M 



2,-2 



2A-3e 



11 

e 



(47r)2 216 
which coincides with Eq. (A. 30) of [2]. 



1^ + 12,. -10^ + 64^ 
2 C(-3) C(-l) 



, (4.21) 
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5 Conclusions 



We have re-examined the three most prominent cases of massless bosonic three-loop vacuum 
sum- integrals, in order to simplify their derivation and translate the original calculations to 
a language that is amenable to generalizations. 

First, we have re-derived the result for the spectacles-type 3-loop vacuum sum-integral 
given first by Andersen and Kyllingstad in [10] , streamlining the computation quite a bit by 
using our notation from [13]. As an improvement over [10], we give a one- dimensional integral 
representation of A (which was given as a triple integral there). Further effort would be 
welcome in order to derive a high-precision result for the numerical coefficient C, involving an 
infinite sum and a one-dimensional integral, leading to extremely slow convergence behavior. 
It would be interesting to study generalizations of the computation outlined in Sec. 2, such as 
1/P^ — 7> 1/[P^]^ as was done for the 3-loop basketball topology in [13], or including factors of 
Pq or other scalar products in the numerator, in order to derive some of the integrals needed 
in our 3-loop computations. 

Second, we have re-derived the result for the basic basketball-type 3-loop vacuum sum- 
integral given first by Arnold and Zhai in [8], streamlining the computation quite a bit by 
using our notation from [13]. 

Third, we have re-derived the result for the first non-trivial 3-loop vacuum sum-integral 
involving scalar products in the numerator given first by Arnold and Zhai in [8], somewhat 
streamlining the computation. Here is a summary of this computation in a nutshell: 

M2,-2-mf = Q,Q.R,R. ^''^^ i^,.? [noo,n^^]2 [U,,]^ + 0{e)r^Mo,o + 0{e) 



One wonders whether there is a simpler way to compute 7W2,-2- Note that the projection 
method, acting on the level of sub-integrals, seems to over-complicate the computation by 
triggering factors of 1/p^ (stemming from = l/P^p'^), which leads outside the class of 

integrals Eq. (1.1) we started with. It even leads outside the natural generalization of this 
class as suggested by IBP methods (which allows for factors of qq etc in the numerators [5]). 
One idea to avoid this change of structure could be to explore applicability of the generic 
tensor method of Ref. [17] to the case of finite-temperature sum-integrals as discussed here. 
However, this is clearly beyond the scope of the present paper but should be explored in the 
future. 

In closing, we hope that our unified exposure of known techniques for sum-integral eval- 
uation leads to a program of generalizing them to other cases - be it with irreducibles in the 
numerator or with different powers of the denominators - as needed for example for deter- 
mining matching coefficients in effective field theories [11], or for advancing to the next loop 
order [9]. In the short term, it seems that the class of hot bosonic tadpoles AiN,-2 is a suit- 
able candidate deserving further study. Finally, an extension to fermionic cases (ultimately 
involving masses as well as chemical potentials) would be another possible line of future work. 
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A Standard integrals 

For convenience, we collect here the functions used above, as defined in [13]. They are the 
1-loop massless propagator at zero temperature 

the 1-loop bosonic tadpoles 



ra_^ Igor _ 2TC{2s-a-d) r{s-l) .^^^rO o^ 



and a specific 2-loop tadpole 

A{si,S2,S3) = A{si,S2,S3;0) (A.3) 



(A.4) 



PQ [Q2]^1[P2]..[(P + Q)2]S3 
_ 2T^ C(2S123 -2d- 54) r(si3 - f )r(si2 " f )r(f - Si)T{si23 - d) 

(27rr)2«i23"2d-s4 (47r)'^r(s2)r(s3)r((i/2)r(sn23 - d) 

where Sahc... = Sc + Sh + Sc + ■■■ ■ 

B Numerical evaluation of A, B, C 

The integrals A and B are easily evaluated numerically e.g. with Mathematica[19], 

A « 9.5763057898113125 , (B.l) 
B ^ 0.058739245719225247 , (B.2) 



while C is tougher to get with high precision. Maybe it is easier to handle it in pieces: 



Ei(-2nx) + In ( — ) 



+0.016232689597 , 



C72 = -6 —( coth(x) - - - V [Ei(-2nx) - e'^"^ In (2nx e^^] 
Jo X \ X 3J f^^^ 



n=l 

-0.022965150204 , 
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C3 = -6 [fcoth(rE) - i 



J Q X \^ X 




1 



+ — ln(4x2e^E 



90 



« -0.021888498587 , 



Jo a; 90 ^ 
« +0.032117000867 , 

C7 = Ci + ^2 + Cs + ^4 



7E + 21n(2)-l 
30 



+0.003496041673 . 



(B.3) 



The sum in Ci converges reasonably fast, the one in C2 rather slowly. The following piece of 
Mathematica code was used to obtain the approximate numerical results given above: 

wp = 60; ag = 30; max = 10000; 

(*liere the sum converges reasonably fast*) 
Clear [ccl] ; 
ccl[n_] : = 

ccl[n] = NIntegrate [-6/ 

X (Coth[x] - 1/x - x/3) (ExpIntegralEi [-2 n x] + 
Exp [-2 n x] Log [2 x/n] ) , {x, 1, Infinity}, 
WorkingPrecision -> wp, AccuracyGoal -> ag] ; 
cl = Sum[ccl[n], {n, 1, max}] 

(*this is the part in which the sum converges extremely slowly*) 
Clear [cc2] ; 
cc2[n_] : = 

cc2[n] = NIntegrate [-6/ 

X (Coth[x] - 1/x - x/3) (ExpIntegralEi [-2 n x] - 
Exp [-2 n x] Log [2 n x Exp [EulerGamma] ] ) , {x, 0, 1}, 
WorkingPrecision -> wp, AccuracyGoal -> ag] ; 
c2 = Sum[cc2[n], {n, 1, max}] 

(*these are not problematic at all*) 
c3 = 
NIntegrate [-6/ 

X ((Coth[x] - 1/x - x/3) /(Exp [2 x] - 1) + x-2/90) (2 Log [2 x] + 
EulerGamma), {x, 0, 1}, WorkingPrecision -> wp, 
AccuracyGoal -> ag] 
c4 = Integrate [x/15 (2 Log[2 x] + EulerGamma), {x, 0, 1}] 

(*sum up to get C*) 



- 14 - 



-Ccl, c2, c3, N[c4, wp]> 
Total [y.] 



C Integrals needed for D, E, F 

Following Appendix C of Ref. [8], let us introduce - in a manner similar to dimensional 
regularization - a convergence factor into convergent integrals of the form Eqs. (3.8), 
(3.10) and (3.13), enabling us to integrate term by term and letting 5 in the end. The 
basic relations needed are 



/>oo 

H{a,h) = / dxx^^^ coth*(2;) 
Jo 

H{a,0) = 



^fa,l) 



H{a,h) 



dx x 



dx x' 



a+5 



i+S 



1 + 2J2' 



-2nx 



n=l 



r(l + a + 5)C(l + a + 5) 



2a+5 



a + 5 
6-1 



H{a-l,b-l) + H{a,b-2) 



(C.l) 
(C.2) 

(C.3) 

(C.4) 



allowing for a recursive solution of integrals Jq°° dx x^ coth^(x) for 6 € Z. 



D Derivation of Eq. (4.14) 

In order to derive Eq. (4.14), we closely follow Eqs. (H. 15-27) of Ref. [8]. The (inverse) 3d 
spatial Fourier transforms that we need derive from Eq. (22) of [13], 



r{s) T fd^ 



(4,) 



E ^"?; (D,i) 



2* (27rr)2* ^ (s-nV. 2"n! 

^ ' n=0 



[P^Y (4^)^ J 

where s = |3/2 - s\ - 1/2, and read (using /3/r(e) = l/i^nf + 0(e) and h = + 0(e)) 

Jpr „-\po\r^ 



2f f 

r2' 3 



1 



d'^r 



e e 



X <! ^, ^, IpoI + ^, -2(2^r)2 + ^ + 1^ ) , e-(l^ol+ko+Poh|Pol)r + 0^^) , (D_2) 



r r 



From Eqs. (4.16) and (4.6) we can therefore compute the transform of Eq. (4.16): 



+ Po 



Q2 (Q + p)2 



IIb 2A 

d P2 



+ •5^0 4/2 



r /■ d^ 



(47r)2./ r2 



1^ ipr^-bo|r(-27^2")2 



e e 



|Po|^^^+PoCoth(f) + 2 



coth(r) + 



sinh (r) 



+ 



+ 



Po , ^IpoI _3 

_ ~r _o ~r -Q 



+pI 



r2 

5,,— + 0(e) 
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{Airy 



^e^P'-a^e-IPol^' ( coth(r) 



7 / / 

-—J dre^lP"!'" ( coth(r 



1 r 
f 3 
f 



rp2 

+ 0(e) 



sm(pr) 



(D.3) 



where in the first hne we have transformed the numerator of the sum-integral in Eq. (4.6) 
as (2go + Po)^ = '^Qo + Po + 2(^0 + Qo)Po — ^ '^Qo — Vq (shifting half of the term linear in go as 
Q — > —Q — P using the denominator's symmetry) and re-arranged terms; in the second line 
we have applied Eq. (D.2), solved the Matsubara sum^^ and used /| = — -|- 0{e); as a 
third step we have done an identical re-writing in terms of derivatives; and in the last line 
we have integrated over angles via / ^ e'P"" = Att J^dr ^^^j^^^ and integrated by parts (twice; 
surface terms cancel against Sp^T"^ /Q). Analogously, for Eq. (4.15) one easily gets 

T 



(47r)2 



-.-{l-5p,)-]+0{e) 



— / dre~lP°l'' ( coth(f) 



r 6 J pr 



(D.4) 
(D.5) 



To prove Eq. (4.14), let us now use Eq. (D.4) on its rhs and integrate over angles: 



Ha 



r2 



PO 



(4vr) 



d^i 



-2\po\r 



T 

Att 



E 



oo 
dr 



■|po|»'. 



Airr 



coth(r) 
coth(f) 



1 
f 

1 

f 



(l-'Jpo) 



+ 0{e) 



+ 0{e) 



(D.6) 
(D.7) 



On the other hand, plugging Eqs. (D.3) and (D.5) into the Ihs of Eq. (4.14) results in 



1 ^ 



d-2 ^ (47r 



dr/dse-|P» I (coth(r 
h 



'l-'^Po)Jj (coth(s) 



(1-5, 



d'^p 1 q2 sin(pr) 



(27r)3 p' 



3 r.2 



■pr 



4 ^' + 2 



sin(j>s) 



ps 



3 

(D 



Using Eqs. (H.21),(H.22) of [8], the integral in the last line of Eq. (D.8) is 

(5(r — s) 



1 
2^ 



d e(r 

+ 2^ 



3rl 



+ {d-2) 



Anrs 



1 1 



2tt d 3rl 



Anrs 



(D. 



where r> = r 6{r—s)+s 9{s—r) and in the last step we used that the first line of Eq. (D.8) is 
symmetric under r o s. Up to terms of 0{e), the second line of Eq. (D.8) thus reduces^^ to 
a delta function, transforming Eq. (D.8) into Eq. (D.7), which completes the proof. 



"From Eq. (19) of [13], it immediately follows that J^qg e"<l''"l + '«'+''«l"lf"'''' = [coth(f) + and 
that 2g^e-(l™l + l™+™l-lf"l)'- = [|po|(l + 2pg)/3 + pgcoth(f) + (coth(f) + \po\) / sinh^ (f)] . 
^^This is the "amazing cancellation" the authors of Ref. [8] refer to in their Appendix H. 
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